Constraints on f(R) Cosmology in the Palatini Formalism 
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In this work we derive the covariant and gauge invariant perturbation equations in general theories 
of f(R) gravity in the Palatini formalism to linear order and calculate the cosmic microwave back- 
ground (CMB) and matter power spectra for an extensively discussed model, f(R) = R + a(— R) 13 , 
which is a possible candidate for the late-time cosmic accelerating expansion found recently. These 
spectra are discussed and found to be sensitively dependent on the value of /3. We are thus able 
to make stringent constraints on from cosmological data on CMB and matter power spectra: 
The three-year Wilkinson Microwave Anisotropy Probe (WMAP) data alone gives a constraint 
\f3\ < 0(1O~ 3 ) while the joint WMAP, Supernova Lagacy Survey (SNLS) and Sloan Digital Sky 
Survey (SDSS) data sets tightens this to j3 ~ 0(1O -6 ), about an order of magnitude more stringent 
than the constraint from SDSS data alone, which makes this model practically indistinguishable 
from the standard ACDM paradigm. 

PACS numbers: 04.50.+h, 98.70.Vc, 98.65.-r 



I. INTRODUCTION 



It is observed that the universe is now undergoing ac- 
celerating expansion [H, 0, H|, which is also consistent 
with the three-year Wilkinson Microwave Anisotropy 
Probe (WMAP) data [3| and several other cosmological 
observations. The usual "explanation" for this involves 
a mysterious component, called the dark energy, which 
drives this accelerating expansion. However, this dark 
energy problem could also be attacked by modifying the 
theory of gravity so that it departs from the standard 
general relativity (GR) when the spacetime curvature be- 
comes small. In one type of modified gravity theories, the 
Ricci scalar R in the Einstein-Hilbert action is simply re- 
placed by a function of R, commonly known as f(R). 
Indeed, in it was shown that by adding correction 

terms, such as R 2 , R ab R a b and R abcd R abcd , to the action, 
the late time accelerating cosmic expansion could be re- 
produced (see also Q and references therein for related 
works). Another argument in favor of such generaliza- 
tions is that the effective Lagrangian of the gravitational 
field generally will include higher order terms in the cur- 
vature invariants as a result of quantum corrections (see, 
e.g., 1). 

However, the conventional metric approach to f(R) 
gravity leads to fourth order equations which may exhibit 
violent instabilities when matter is present in the weak 
gravity regime 0] (see however [l(| for a discussion) . On 
the other hand, in the Palatini variation of the action 
where the metric and connection are treated as indepen- 
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dent dynamical variables [ll| , the resultant equations are 
second order, which are more tractable and concordant 
with field equations in other branches of physics. In par- 
ticular, the typical form f(R) — R + a(-R) 13 has been 
discussed extensively in the literature as an alternative 
dark energy model which fits rather well with the su- 
pernovae (SNe) la data, and it is also tested using cos- 
mic microwave background (CMB) shift parameter and 
baryon acoustic oscillation (BAO) in [l2|, 0, EH [H| ■ Pos- 
sible constraints on other types of Palatini-/ (R) model 
have also been considered using big bang nucleosynthesis 
(BBN) and the requirements that the success of the infla- 
tionary paradigm is not spoiled [l6| ■ Recently constraint 
from data on matter power spectrum alone is given in 
[ItJ • As far as we know, there have been no attempts to 
confront Palatini f(R) gravity models with CMB data to 
date. 

In this work, we will concentrate on the model of 
f(R) = R + a(—R)P, where a is positive (so that it can 
reproduce the recent cosmic acceleration). We refer to it 
as the late f(R) cosmological model because its correc- 
tions to GR dominate very lately, and we study both its 
CMB and matter power spectra. For this, we need the 
perturbation equations in the Palatini formalism, one set 
of which has been worked out in (l8j . Here, however, we 
will derive a set of covariant and gauge invariant pertur- 
bation equations by the method of 3 + 1 decomposition 
(see Sec. [TTJ) for our calculations. Also we shall try to 
constrain the model parameters. Unlike previous works, 
we use the full three-year WMAP data set [| instead of 
the CMB shift parameter only. This f{R) gravity model 
will be constrained firstly by the WMAP CMB spectra 
data, and then jointly by the CMB spectra, SNe mea- 
surements from Supernova Legacy Survey (SNLS), plus 
the matter power spectrum data measured by the Sloan 
Digital Sky Survey (SDSS) [fl. 
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The paper is organized as follows. In Sec. HH we will 
first review briefly the theory of f(R) gravity in the Pala- 
tini formalism, and we then present the perturbation 
equations. Then in Sec. IHII the CMB and matter power 
spectra for different choices of parameters are displayed 
and discussed, and the constraints from various data sets 
will be presented. Finally, we conclude in Sec. lIVI There 
is also an appendix where we show that our perturba- 
tion equations are equivalent to those in the synchronous 
gauge under specified conditions. Throughout this work 
we will assume a flat universe filled with cold dark mat- 
ter, photons, baryons, electrons and 3 species of neutri- 
nos (all massless); the unit c = 1 is adopted. The metric 
convention used in this paper is (+,—,—,—). 



Eq. © reads 



FR — 2/ = kT 



(4) 



with T = p — 3p (p is the energy density and p 
the isotropic pressure) being the trace of the energy- 
momentum tensor. This is the so-called structural equa- 
tion [2l[ which relates R directly to the energy compo- 
nents in the universe: given a specific form of f(R) and 
thus F(R), R can be obtained as a function of T by nu- 
merically or analytically solving this equation. 

The variation of Eq. ([T]) with respect to the connection 
field f leads to another equation 



V a [F{R)^—gg bc ] = 0, 



(5) 



II. FIELD EQUATIONS IN THEORIES OF 
PALATINI - f(R) GRAVITY 

In this section we shall first summarize the properties 
of the general theory of f(R) gravity in the Palatini for- 
malism, and we then derive its perturbation equations. 
These equations could be found elsewhere [2(| , and here 
we list them just for completeness. 



A. General theory of f(R) gravity in the Palatini 
approach 

The starting point of our discussion on the Palatini- 
f(R) gravity is the modified Einstein-Hilbert action, 
which is given as 
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f(R) + Cn 



(1) 



where k = 8irG (G is the Newton's constant) and R = 
g ab R ab (f) (a, 6 = 0,1, 2, 3) with R ab (T) being defined as 
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(2) 



Notice that the connection T here is not the conventional 
Levi-Civita connection of the metric g ab , which we shall 
denote by T; rather it will be treated as an independent 
field in the Palatini approach to the f(R) theory of grav- 
ity. Correspondingly, the tensor R ab and scalar R are 
also not the Ricci tensor and Ricci scalar calculated from 
g ab as in GR, which instead are denoted by lZ ab and 1Z 
respectively in this work (1Z — g ab 7Z ab )- The matter La- 
grangian density C m , on the other hand, is assumed to 
be independent of T, which is the same as in GR. 

The extremization of the action Eq. ([T]) with respect to 
the metric g a b then gives the modified Einstein equations 



FR a b 
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(3) 



in which F = F(R) = df{R)/dR and T ab is the energy- 
momentum tensor in the system discussed. The trace of 



which indicates that the connection T is compatible with 
a metric j ab that is conformal to g ab : 



7a h = F(R)g ab 



(6) 



With the aid of Eq. |6j) we could now easily obtain the 
relation between R ab and lZ ab as 
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Note that in above we are using T> and V to denote the 
covariant derivative operators which are compatible with 
g ab and j ab respectively. 

Since C m depends only on g ab (and, of course, some 
matter fields) and the energy-momentum conservation 
law holds with respect to it, we shall treat this metric 
as the physical one. Consequently the difference between 
f(R) gravity and GR could be understood as a change in 
the manner in which the spacetime curvature and thus 
the physical Ricci tensor lZ ab responds to the distribution 
of matter (through the modified Einstein equations) . In 
order to make this point explicit, we can rewrite Eq. ([3]) 
by the use of Eq. ([7]) as 

K% b = FTZ ab - -g ab f 

+ ^V a FV b F - V a V b F - X -g ab V c V c F, (8) 
in which F(T), f{T) are now simply functions of T. 



B. The Perturbation Equations 

The perturbation equations in general theories of f(R) 
gravity have been derived in [18j ]. However, here we 
adopt a different, covariant and gauge invariant deriva- 
tion which utilizes the method of 3 + 1 decomposition 

The main idea of 3+ 1 decomposition is to make space- 
time splits of physical quantities with respect to the 4- 
velocity u a of an observer. A projection tensor h ab is 
then defined as h ab = g ab — u a u b which could be used to 
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obtain covariant tensors orthogonal to u. For example, 
the covariant spatial derivative D of an arbitrary tensor 
field T^.'.'g (which, by definition, is orthogonal to u) is 
given as 



V a T b-c 



■ Kh r d - 



hrn'Tc; 



k 

i — s ' 



(9) 



The energy-momentum tensor and covariant derivative 
of the 4-velocity u could be decomposed respectively as 



Tab = K ab + 2q {a u b) + pu a u b - ph ab , 



V a u b 



Oab 



ZVab + ~0h ab + U a A b . 



(10) 

(11) 



In the above ir ab is the projected symmetric trace free 
(PSTF) anisotropic stress, q the vector heat flux, a ab the 
PSTF shear tensor, w ab = T>\ a u b ], 9 = T> a u a = 3d/ a 
(a is the cosmic scale factor) the expansion scalar and 
A a = u a the acceleration. The overdot expressed as 
(j> = u a T> a <j) is the derivative with respect to the proper 
time of the comoving observer moving at velocity u, 
and the square brackets denote antisymmetrization and 
parentheses symmetrization. The normalization is cho- 
sen to be u 2 = 1 in consistence with our metric conven- 
tion. 

Decomposing the Riemann tensor and making use of 
the modified Einstein equations with the general tech- 
niques used in GR, we obtain, after linearization, five 
constraint equations 
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Here e abc d is the covariant permutation tensor, £ ab and 
B a b are respectively the electric and magnetic parts of 
the Weyl tensor W a bcd, given respectively by £ ab = 
u c u d Wacbd and B ab = ~±u c u d e ac ef W efbd . 
We also obtain seven propagation equations: 

P + {p + P)0 + t> a q a 



4a + ^0q a + (p + p)A a - VaV + V Kab 



0;(17) 
0;(18) 
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The angle brackets mean taking the trace free part of a 
quantity. 

Besides the above equations, it would also be useful to 
express the projected Ricci scalar 1Z in the hypersurfaces 
orthogonal to u a (the projected Riemann tensor, lZ a bcd, is 
defined by [V a , V b \v c = TZ abd c v d , similar to the definition 
of the full covariant Riemann tensor TZ abcd but with a 
conventional opposite sign, and the calculations for the 
projected Ricci tensor !Z a b and projected Ricci scalar 1Z 
just follow the same way as in GR) as 
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k(p + 3p) - / 
F 



3F 
2F 



2V 2 F 
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The spatial derivative of this projected Ricci scalar, rj a = 
^aD a Ti-: is then obtained after a lengthy calculation as 

" 3 
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(25) 



and its time evolution is governed by the propagation 
equation 

26» 



Va + Y Va 
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2F 
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F 
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As we are considering a spatially flat universe, its spa- 
tial curvature will vanish for large scales, meaning that 
It = 0. Thus from Eq. flU} we have 



F 
2F 
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6F 



[ K (p + 3p)-f]. 



(27) 
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This is just the modified (first) Friedmann equation in 
the f(R) version of gravitational theory, and the other 
modified background equations (the second Friedmann 
equation and the energy-conservation equation) could be 
obtained by taking the zeroth-order parts of Eqs. (fl7|) 
and (fT9|) . It is easy to check that when f(R) = R, we 
have F = 1, and these equations just reduce to those 
in GR - in this case GR and the Palatini-/(i?) theory 
lead to the same results. In the appendix we also show 
that this set of perturbation equations is equivalent to 
that derived in the synchronous gauge, which serves as a 
check for this work. 

Remember that we have had /, F and R as functions 
of T at hand, it is then straightforward to calculate 
F,F,V a F,V a F etc. as functions of f = -(p b + p c )6 
and T> a T = (1 — 3c 2 )T> a pi, + V a p c , in which p b ^ is the 
energy density of baryons (cold dark matter) and c s is the 
baryon sound speed. Note that in this work we choose 
to neglect the small baryon pressure except in the terms 
where its spatial derivative is involved, in which case they 
might be significant at small scales. The above equations 
could then be numerically propagated given the initial 
conditions, to obtain the evolutions of small density per- 
turbations and the CMB and matter power spectra in 
theories of f(R) gravity. Finally the three- year WMAP 
data on CMB spectra, SNLS SN data and SDSS data on 
matter power spectrum could be used to constrain pa- 
rameters in the f(R) models. These results will be given 
in the following section. 



III. NUMERICAL RESULTS AND 
COSMOLOGICAL CONSTRAINTS ON THE 
MODEL 

This section is devoted to numerical results and con- 
straints of the present model. To this effect we will 
first very briefly summarize and explain the effects of 
the f(R) modifications to GR on the linear spectra; for 
more details see [2(| ■ After that we shall employ the pub- 
lic Markov Chain Monte Carlo (MCMC) engine [H] to 
search the parameter space with the theoretical CMB and 
matter power spectra calculated by the modified CAMB 
code; the constraints are then summarized and discussed. 

In Fig. Q] we have displayed the TT CMB spectra for 
the model with different choices of (3. It is obvious from 
this figure that, when f3 < 0, the spectrum gets a boost 
in the scales I < 100, which could be significant if \(3\ 
is large enough. This effect is due to a strong late-time 
integrated Sachs- Wolfe (ISW) effect [H HI, which in 
turn originates from the unusually rapid late-time decay 
of the gravitatoinal potential cf> of the present f(R) model 
compared with ACDM, as shown in the lower panel of 
Fig. 3 in [2(| (see this reference for more details). We 
have not given the curves for (3 > because in that case 
the spectrum generally blows up except for very small 
|/3|s (see below). 

Another interesting feature in Fig.[T]is that for negative 




FIG. 1: The TT CMB spectrum for the f(R) = R + a(-R) 13 
model, with Q. m (current fractional energy density of nonrela- 
tivistic matter) and Ho (current Hubble constant) fixed to be 
0.3 and 72 km/s/Mpc respectively. Choices of (3 are indicated 
besides the curves. The case /3 = corresponds to a ACDM 
Universe. 



(3 the spectrum shifts towards the right-hand-side (larger 
Vs), likely due to the unusual angular-distance-redshift 
relation [2^, [2?| ■ Since the standard ACDM cosmology 
is expected to be valid in the early times when the cor- 
rection to GR is negligible, the sound horizon and the 
thickness of the last scattering surface are the same as 
in ACDM. But at late times the Friedmann equation is 
modified (Eq. (|2~T)) ), and so is the relation between red- 
shift and conformal distance. This would cause the CMB 
spectrum to shift sideways. This shifting effect, however, 
is negligible for the constrained ranges of j3 obtained be- 
low. 

The CMB EE polarization and cross correlation spec- 
tra show no additional interesting features and cannot be 
used to put strong constraint on the model parameters, 
and so we will not present and discuss them here. 

We have also given in Fig. [J] the matter power spectra. 
As indicated in this figure, the matter power spectrum 
depends sensitively on the value of (3 and could differ 
from ACDM significantly even if \f3\ only deviates from 
by a tiny amount e.g., of order 0(1O~ 5 ). This fea- 
ture has been pointed out and discussed extensively in 
El, H3| • Basically, this is because of the sensitive re- 
sponse of the modified gravity to the spatial variations of 
matter distribution, which, at small enough scales would 
significantly affect the growth of density perturbations. 
To be explicit, for large enough k's, the growth equation 
for the comoving energy density fluctuations S m could be 
written as (l7| 



d 2 S m 
dN 2 



F 



a 2 H 2 3F(2FH + F) 



in which N = log(a) and F/3F(2FH + F) = 



%eff 



(28) 



acts 



FIG. 2: The matter power spectra of the f(R) = R+a(-Rf 
model for different choices of /3 (both negative and positive) 
as indicated beside the curves. The case (5 = corresponds 
to a ACDM Universe. 



FIG. 3: The marginal distributions of Q m and (3, obtained us- 
ing the three year WMAP data alone. Here the distributions 
are normalized such that the maximum probability density is 
1. 



as an effective sound speed squared that vanishes in the 
ACDM model. For the present f(R) model, we have 
F = 1 — a(3(— R)"^ 1 , in which a and —R are positive 
and — R decreases with time. So if (3 < 0, then F > 
and thus c 2 e ^j > 0; this effective pressure term will 
restrict the growths of small-scale density perturbations 
and leads to oscillations (as shown in Fig. [3 for the case 
of (3 = —0.00005) of the spectrum in these scales. On 
the other hand, if (3 > (as we hope to recover standard 
ACDM cosmology in earlier times, we shall also restrict 
ourselves to /3 < 1), then F and c 2 e ^ will be negative; 
this will make the density fluctuations unstable and blow 
up, the same reason why the CMB spectrum depends so 
sensitively on positive /3s. 

Since from Figs. [1] and [5] we have seen that the linear 
spectra of our f(R) model depend very sensitively on 
the model parameter (3, it can be expected that the data 
on CMB and matter power spectra could place stringent 
constraints on (3, as we will show now. As mentioned 
in Sec. HI we shall firstly use the full three year WMAP 
data and then perform a joint constraint simultaneously 
using WMAP, SNLS and SDSS data to constrain the 
parameters. 

Because the Hubble parameter Ho is already measured 
to rather good precision by the HST Key Project, we 
shall use Hq = 72 km/s/Mpc in our calculations. 
Therefore we vary the following parameters: baryon den- 
sity uif, — flbh 2 , cold dark matter uj c — fl c h 2 , reioniza- 
tion rcdshift z re , spectral index n s , normalization am- 
plitude A s and the model parameter (3. In Fig. [3] the 
marginal distributions of £l m and (3 are shown. The 95% 
confidence interval for £l m and (3 are [0.233, 0.268] and 
[—3.45 x 10~ 3 , 3.07 x 10~ 3 ] respectively. We also present 
the contour plot of the joint distribution of tt m and (3 
in Fig. [3J From these figures we can see that the CMB 



at 0.25 




0.5 

0.45 

0.4 

0.35 

0.3 

0.25 

0.2 

0.15 

0.1 

0.05 





FIG. 4: (Color Online) The contour plot of the joint distri- 
bution of Q m and j3, constrained by WMAP data alone. The 
inner and outer loops are the 68% and 95% confidence con- 
tours respectively. 



spectra could constrain \(3\ to 0(1O -3 ), ~ 100 times more 
stringent than the constraint from the CMB shift param- 
eter [lj], which is of order 0.1. That the CMB spectra is 
much more powerful in constraining the parameters than 
the CMB shift parameter is expected because the former 
bears a lot more information than the latter. It looks 
from these figures that a slightly positive (3 is preferred 
by the CMB data. 

To tighten the bounds on the parameters, we per- 
form a joint constraint making use of both the three- 
year WMAP and the SDSS data sets. In addition to the 
matter power spectrum, the SNe data from SNLS [3(| is 
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FIG. 5: The marginal distributions of the various model parameters, constrained simutaneously by the WMAP, SNLS and 
SDSS data sets. The distributions are normalized such that the maximum probability density is 1. 

also used in the joint constraint, though their effects are 
found to be negligible. For the SDSS data, we conser- 
vatively adopt the measurements for scales larger than 
k = 0.2h Mpc -1 (where h = H /(100 km s^Mpc" 1 )) to 
avoid encountering the nonlinear effects in the measured 
matter power spectrum. The bias between galaxy power 
spectrum and matter power spectrum is assumed to be 
a scale independent constant; CosmoMC [28| assumes a 
flat prior on it and marginalizes analytically. 

The calculation indicates that indeed the allowed range 
is shrunk, as indicated in Figs. [5] and [5] (for completeness 
in Fig. [7] we have also plotted the best-fitted curves with 
the observational data points from SNLS, WMAP and 
SDSS we use in the constraints). The 95% confidence 
intervals for Q m and (3 now become [0.241, 0.274] and 
[2.12 x 10~ 6 , 5.98 x 10" 6 ] respectively, and the 95% con- 
fidence interval for a% is [0.85, 1.21]. The distribution of 
Q m does not change much since it is already well con- 
strained by the WMAP, but the bound on (5 is tightened 
to the order of 10~ 6 (and obviously future refined data 
could still further this constraint). What is more, these 
joint constraints also prefer a positive (3 and actually have 
excluded the case of f3 = 0, i.e., the ACDM paradigm, 




FIG. 6: (Color Online) The contour plot of joint distribution 
of rim and /3 under the constraints of the WMAP, SNLS and 
SDSS data sets. The inner and outer loops are the 68% and 
95% confidence contours respectively. 
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FIG. 7: (Color Online) The data points from SNLS (upper 
panel), WMAP (middle panel) and SDSS (lower panel) data 
sets against the theoretical curves of our best-fitted model. 
For SDSS data we have plotted both vertical and horizonal 
error bars. Note that the last three data points from SDSS 
(for which k > 0.2h Mpc -1 ) are not included in our numerical 
constraints, and we have used a bias of 1.1 to relate the best 
fit theoretical power spectrum to the SDSS data. 

at the 95% confidence level. In Fig. [7] we could see that 
all the data sets are fitted very well. Furthermore, more 
stringent constraints on [3 can be obtained if data points 
in the nonlinear regime are also used since the matter 
power spectrum will blow up in the small scales for posi- 
tive (3 (see Fig. Nonetheless, our stringent constraint 
on (3 already makes the model nearly indistinguishable 
from ACDM for practical purposes, and without a natu- 
ral motivation for such tiny values of (3 this model should 
be more reasonably disfavored. 



IV. DISCUSSION AND CONCLUSION 

In conclusion, we have in this work derived the per- 
turbation equations for general theories of Palatini-/^) 
gravity and applied them to a typical class of model 
f(R) = R + a(—R)° which is proposed as an alternative 
to the cosmological constant to account for the late-time 
accelerating cosmic expansion and has been extensively 
studied. We then calculate the CMB and matter power 
spectra for this model using a modified CAMB code. It 
is shown that for negative /3s the potential <fi will see an 
unusually rapid decay at late times, leading to an en- 
hancement of the ISW effect and thus a boost of the TT 
CMB spectrum at small Is. There also appears a posi- 
tive effective pressure term in the equation governing the 
growth of density perturbations, which could be signifi- 
cant for small scales (large ks) and restricts the pertur- 
bation growths in these scales. For positive /3s, however, 
the small-scale density fluctuations will become unstable 
and grow exponentially, resulting in blowing-ups of the 
matter power spectrum. 

We have constrained the model parameters using the 
WMAP, SNLS and SDSS data. Because the CMB and 
matter power spectra are rather sensitive to the exact val- 
ues of the parameter /3, we are able to give much more 
stringent constraints on /3 (0(1O -3 ) and 0(1CP 6 ) respec- 
tively) than those (©(lCT 1 )) coming from the CMB shift 
parameter fitting or measurements on SNIa [l3[. Com- 
pared with the bound (0(1CT 5 )) from SDSS data alone 
our WMAP + SNLS + SDSS constraint is tighter 
because the allowed range of fl m is largely reduced here. 

These constraints seem to make the present model (in 
its allowed parameter space) indistinguishable from the 
ACDM paradigm and raise a fine-tuning problem to the 
late f(R) gravity theory. However, there still remains the 
interesting possibility that f(R) modification of gravity 
enters at earlier times (or higher densities): can it survive 
the tests from WMAP and SDSS data? This topic is 
beyond the scope of this article and has been investigated 
in another work [20] ; in that case the parameter space for 
f(R) gravity is also highly limited. 

Appendix 

This appendix is devoted to showing the equivalence 
between the covariant perturbation equations derived 
here and the perturbation equations in the synchronous 
gauge [3l| . In the latter case, the line element is ex- 
pressed as (indices i, j, k run over 1, 2, 3 hereafter) 

ds 2 = a(r) 2 [dr 2 - (fly + hf^dx'dx^ , (29) 

in which r is the conformal time given by dt = a{r)dr 
and we have defined the scalar modes in Fourier space as 
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4(x,r) 



/ 



d 3 k exp(ik • x) 



kikjh s (k,T) + 6 



(30) 



where a superscript 'S" denotes quantities in the syn- 
chronous gauge and k = k/fc is the unit vector in the 
k-direction. 

In order to relate the synchronous gauge variables 
h s and r/ s with those in the covariant perturbation 
method, let us do the harmonic expansion for the first- 
order quantities h a = X> a a and r\ a as h a = J^k khQ k 
and Va = J2 k ^Q k (Here Q k = fV a Q k and Q k is 
the eigenfunction of the generalized Helmholtz equation 
a 2 V 2 Q k = k 2 Q k . For more details see e.g., @, 111), 



to obtain the variables h and r\ relevant to specified k- 
modes. We shall choose the reference velocity u to be 
the 4-velocity of cold dark matter, which means that the 
cold dark matter heat flux q c = p c v c and the acceler- 
ation A both vanish [2^, f2^| . Then using the relations 
rj S = —rj/2 and h' s = 6h' [32| (a prime here means tak- 
ing derivative with respect to the conformal time r) it 
is easy to show the equivalence between these two sets 
of perturbation equations. More explicitly: Eq. ([26]) is 
equal to the following first-order equation 



Fk 2 V ' s = ^a 2 (p s +p s )9 s 



3F' 
~2F 



k 2 SF + -k 2 SF', 



(31) 



where TL = a' /a. Taking the spatial derivative of Eq. |T 
leads to 



F 
~2 



h" s + Hh' s + ^h' s 
2F 



2 a 2 df(f) xrr 3F' 36F" 

na dp H —^dl =rdF H 

H 2 dT F 2 



3a'' 



3 / F' 



2 V F 



61(32) 



Then from Eqs. (|25|) and ([32| we can obtain 



F 



,, 2 5 5Hh' s h" s 5F' /S 

Ak 2 rf =-h b 

'2 2 4F 



— Qrj^dT - na STi + 6H6F + 



3a' 1 



3H Z 



5k 2 



$F[33) 



and finally, taking time derivative of Eq. ()26|) and making 
use of Eqs. (HU), (J25j>, ^} and (JSU), we get 



J 



F 



h 2 v -\n{h' s 



677 



- (h" s + 6r," s ) 
3 



na 2 (p + p)cr s + \k 2 8F + ~F' (ti s + 6if) 



(34) 



In the above S means the spatial variation of a quan- background value (up to first order all the barred quan- 
tity for the discussed length scale (or k) and a bar its 
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tities in these equations could be equally replaced by 
their unbarred correspondences). The variables 9 ,a 
are defined respectively as [3l[ (p + p)9 s = ik^STj and 
(p+p)a s = -(^^-^^(Tj-i^T^). T ab is the energy- 
momentum tensor in the synchronous gauge (to be dis- 
tinguished from the T a b used above) and T its trace. 

We have checked that Eqs. (JSTJ) - (ED) are just the 
perturbation equations one has in the synchronous gauge 
for theories of Palatini-/ (R) gravity, as expected [23, [2J]. 



cal calculations. The work is partially supported by a 
grant from the Research Grants Council of the Hong 
Kong Special Administrative Region, China (Project 
No. 400707). B. L. acknowledges supports from the Over- 
seas Research Studentship, Cambridge Overseas Trust 
and the Department of Applied Mathematics and Theo- 
retical Physics at the University of Cambridge. We are 
grateful to Dr. Tomi Koivisto for his helpful discussions. 



Acknowledgments 

We would like to thank the ITSC of the Chinese Uni- 
versity of Hong Kong for using their clusters for numeri- 



A. G. Riess et al , Astron. J 116, 1009 (1998). [arXiv: [17 
|astro-ph/9805201| ] 
S. Perlmutter et a ., Astrophys. J 517, 565 (1999). [arXiv: [18 
|astro-ph/9812133| ] 

A. G. Riess et al Astron. J 607, 665 (2004). [arXiv: [19 
|astro-ph/0402512| ] 

G. Hinshaw et al , |as"tro-ph/0603451| G. Hinshaw et al, [20 
|astro-ph/0603450| 

S. M. Carroll, A. De Felice, V. Duvvuri, D. A. Easson, [21 
M. Trodden an d M. S. Turner, Ph ys. Rev. D 71, 063513 
(2005). [arXiv: astro-ph/0410031| ] 

S. Capozziello, V. F. Cardo ne and A. Tr oisi, Phys. Rev. D [22 

71, 043503 (2005). [arXiv: |astro-ph/0501426|] 

S. Nojiri and S. D. Odintsov, arXiv: |hep-th /0601213 [23 

N. D. Birrell and P. C. W. Davis, Quantum Fields 

in Curved Spacetime, Cambridge University Press, UK, [24 

1982. 

A. D. Dolgov and M. Kawasaki Phys. Lett. B 573, 1 
(2003). [arXiv: astro-ph/0307285| ] [25 
S. Nojiri and S . D. Odintsov, P hys. Rev. D 68, 123512 
(2003). [arXiv: |hep-th/0307288| ] 

D. N. Vollick, Phy s. Rev. D 68, 063510 (2003). [arXiv: [26 
|astro-ph/0 306630 .] 

Capozziello, V. F. Cardone and M. Francav- [27 
iglia, Gen. Rel . Grav. 38, 711 (2006). [arXiv: 
|astro-ph/0410135| ] [28 
M. Amarzguioui, O. Elgaroy, D. F. Mota and T. Mul- 
tamaki, Astron. Astrophys. 454, 707 (2006). [arXiv: [29 
|astro-ph/0510519| ] 
S. Fay, R. Tavakol an d S. Tsujikawa, Ph ys. Rev. D 75, [30 
063509 (2007). [arXiv: |astro-pn/0701479| ] 

A. Borowiec, W. Godlowski and M. Szyd- [31 
lowski, Phys. R ev. D 74, 043502 (2006). [arXiv: 
[astro -ph/0602526|] [32 
T. P. Sotiriou, Class Quant. Grav. 23, 1253 (2006). 
[arXiv: |gr-qc/0512017[ ] 



T. Koivisto, Phys . Rev. D 73, 083517 (2006). [arXiv: 
|astro-ph/0 602031 .] 

T. Koivisto and H. Kurki-Suonio, Class. Quant. Grav. 

23, 2355 (2006). [arXiv: |astro-ph/0509422| ] 

M. Tegmark et al , Astrophys. J 606, 702 (2004). [arXiv: 

|astro-ph/03i0725j ] 

B. Li a nd M. -C. Chu, Ph ys. Rev. D, 74, 104010 (2006). 
[arXiv: |astro-ph/06 10486 .] 

G. Allemandi, ~K. Borowiec and M. Francav- 
iglia, Phys. R ev. D 70, 043524 (2004). [arXiv: 
|astro-ph/0 403264 ,] 

G. F. R. Ellis and M. Bruni, Phys. Rev. D 40, 1804 
(1989). 

A. Challinor an d A. Lasenb y, Astrophys. J. 513, 1 (1999). 
[arXiv: |astro-ph/9804301| ] 

A. Lewis 2000, PhD thesis, Astrophysics 

Group & Cavendish Laboratory, University. 

[http://www.mr ao.cam.ac.uk/'^amll005/cmb| ] 

T. Koivisto, H. Kurki-Suonio and F. Ravn- 

dal, Phys. Rev . D 71, 064027 (2005). [arXiv: 

|astro-ph/0409163[ ] 

A. Riazuelo an d J. -P. Uzan, Phy s. Rev. D 66, 023525 
(2002). [arXiv: |astro-ph/0107386| ] 

J. -P. Uzan, Re v. Mod. Phys. 75, 403 (2003). [arXiv: 
|hep-ph/0205"340[ ] 

A. Lew is and~S~Bridle, P hys. Rev. D 66, 103511 (2002). 
[arXiv: |astro-ph/0205436[ ] 

W. L. Freedman et al, Astrophys. J 553, 47 (2001). 
[arXiv: |astro-ph/0012376[ ] 

P. Asti er et al, Astron. Astrophys. 447, 31 (2006). 
[arXiv: |astro-ph/0510447[ ] 

C. -P. Ma and E. Bertschinger, Astrophys. J 455, 7 
(1995). [arXiv: [astr^ph/9401007 ] 

See, A. Lewis 2000, CAM B Notes, for a description, 
http : / /cosmologist .info /notes / ] 



